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ABSTRACT 

Aims. A set of new, sensitive, and spectrally resolved, sub-millimeter line observations are used to probe the warm circumstellar gas around the 
S-type AGB star^ Cyg. The observed lines involve high rotational quantum numbers, which, combined with previously obtained lower-frequency 
data, make it possible to study in detail the chemical and physical properties of, essentially, the entire circumstellar envelope of x Cyg. 
Methods. The data were obtained using the HIFI instrument aboard Herschel, whose high spectral resolution provides valuable information about 
the line profiles. Detailed, non-LTE, radiative transfer modelling, including dust radiative transfer coupled with a dynamical model, has been 
performed to derive the temperature, density, and velocity structure of the circumstellar envelope. 

Results. We report the first detection of circumstellar H2O rotational emission lines in an S-star. Using the high- J CO lines to derive the parameters 
for the circumstellar envelope, we modelled both the ortho- and para-H 2 lines. Our modelling results are consistent with the velocity structure 
expected for a dust-driven wind. The derived total H2O abundance (relative to H2) is (1.1 ± 0.2) x 10~ 5 , much lower than that in O-rich stars. The 
derived ortho-to-para ratio of 2. 1 ±0.6 is close to the high-temperature equilibrium limit, consistent with H 2 being formed in the photosphere. 
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1. Introduction 

Observations of the dust and gas components in the circumstellar 
envelopes (CSEs) around asymptotic giant branch (AGB) stars 
have been carried out at different wavelengths. Observations 
in the infrared trace the dust as well as the warm molecu- 
lar layer close to the s tellar photosphere (e.g.. lJusttanont et al.l 
Il996t lAoki eTai1ll999t ISchoier et al.ll2002h . Submillimeter and 
radio observations of trace molecules h ave been used to study 
the cooler outer part s of th e CSEs (e.g.. | Knapp & Morrislll985t 
ISchoier & Olofssonll2000t iKemper et alj|2003l) . To bridge this 
gap, the Infrared Space Observatory (ISO) was used to ob- 
serve a large number of AGB stars up to almost 200 fim, and 
in O-rich stars, a number of H2O emiss ion lines were detected 
dBarlow et al.lll996t iNeufeld et alj|1996l) . However, the circum- 
stellar far-infrared lines were unresolved. Hence, crucial infor- 
mation about the line profiles remained unknown. 

Water is an important molecule in CSEs as it is thought to be 
one of the main cooling agents in the wind. It is also expected to 
be a good probe of the inner regions of the CSE where the gas 
is accelerated. However, to fully explore the potential of H2O 
lines as a probe of the circumstellar gas, a full radiative transfer 
has to be performed. Owing to difficulties in calculating accurate 
collisional rates coupled with the very high optical depth of the 
H2O lines in the inner region of the CSE, slow progress has been 
made. Nevertheless, calculations of the heating and cooling in 



the CSEs of O-rich stars suggest that H2O dominates the cooling 
in most parts of the envelope until it is p hotodiss ociated by inter- 
stella r UV photons (Goldr eich & Scoville 1976; Justt anont et al.l 
Il994t IMaereker et al.ll2008l l2009b iDecin et alfeoiOal) . and that 
some lines should come mainly from the acceleration zone. 
Eventually, spectrally resolved circumstellar H2O lines were ob- 
served by two space missions dedicated to search for cosmic 
water line emission: SWAS and Odin. Both missions were able 
to detect the grou nd-state line of H2O at 557 GHz in a num- 
ber of AGB sta r s (lHarwit & Berginl 120021: iMelnick et all 1200 lb 
Justtanont et al.1 l2005t lHasegaw a et al.l 120061: IMaereker et al. 



20091) . It was shown that not only the line intensity, but also the 



line profile is crucial for interpreting the data corre ctly. 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



In 2009, the ESA-Herschel Space Observatory dPilbratt et al.l 
l2010l) was launche d with the Heter odyne Instrument for the Far- 
Infrared (HIFU, d e Graauw et al.l 2010), which aims to study 
H2O line emission in different environments in our Galaxy and 
beyond. HIFI offers the opportunity to study the warm molec- 

1 HIFI has been designed and built by a consortium of insti- 
tutes and university departments from across Europe, Canada and the 
United States under the leadership of SRON Netherlands Institute 
for Space Research, Groningen, The Netherlands and with major 
contributions from Germany, France and the US. Consortium mem- 
bers are: Canada: CSA, U.Waterloo; France: CESR, LAB, LERMA, 
IRAM; Germany: KOSMA, MPIfR, MPS; Ireland, NUI Maynooth; 
Italy: ASI, IFSI-INAF, Osservatorio Astrofisico di Arcetri- INAF; 
Netherlands: SRON, TUD; Poland: CAMK, CBK; Spain: Observatorio 
Astronomico Nacional (IGN), Centra de Astrobiologia (CSIC-INTA). 
Sweden: Chalmers University of Technology - MC2, RSS & GARD; 
Onsala Space Observatory; Swedish National Space Board, Stockholm 



K. Justtanont et al.: A HIFI preview of warm molecular gas around^- Cyg 



Table 1. HIFI observations of x Cyg. 



Molecule Transition v (GHz) E™ (K) /"(Kkms" 1 ) 



CO 


7=6-5 


691.473 


116 


15.3 




CO 


7=10-9 


1151.985 


304 


13.6 




CO 


7=16-15 


1841.346 


752 


13.5 


1 1 


o-H 2 


1 10 ~~ loi 


556.936 


61 


3.0 




p-H 2 


2n - 2o2 


752.033 


137 


4.2 


3 


p-H 2 


2o2 _ In 


987.927 


101 


7.3 


E-i 


p-H 2 


In _ Ooo 


1113.343 


53 


8.0 




o-H,0 


3l2 - 2 2 i 


1153.127 


249 


7.5 




o-H 2 


321 — 3 12 


1162.912 


305 


1.5 




o-H 2 


3o3 ~~ 2[2 


1716.770 


197 


18.1 





Note to the table: "The absolute calibration accuracy is between 10% 
to 30% (see text). 



ular layers in CSEs of AGB stars in great detail, e.g., the high 
spectral resolution and wide spectral coverage allow a detailed 
study of the gas dynamics. 

As part of the guaranteed time programme HIFISTARS 
(P.I.: V. Bujarrabal), the S-star (C/O* 1) x Cyg was selected 
for st udy. Distance estimates range from 150 pc dKnapp et al.l 
120031) to 180 pc (Ivan Leeuwenl 12007(1 The star exhibits 



SiO masers (e.g.. lOlofsson et al. 



198lHSchwartzetal]ll982t 



lAlcole a & Buiarra baj 19921). but no H7O maser emi ssion has 



been found ( Mente n & Youngll995l:IShintani et al.ll2008l) . Being 
nearby and bright, x Cyg has bee n observed using in terferomet- 
ric tec hniques in both the opt ical dLacour et al.ll2.Q09b and the in- 
frared dTevousian et alj|2004l) . In this Letter, we briefly describe 
the observations in Sect. 2, we discuss the modelling of the ob- 
served molecular emission in Sect. 3, and present our results and 
conclusions in Sect. 4. 



2. Observations 

The HIFI data were obtaine d using the dual-beam-switching 
mode dRoelfsema et al.ll2.010t> with a throw of 3' and slow (0.5- 
1Hz) chopping in March- April 2010. A total of 8 frequency set- 
tings with a total of 10 lines detected are being reported in this 
paper. The targeted lines were selected to cover a wide range of 
excitation temperature, exploring different regions of the CSE. 
As backend, the wide-band spectrometer (WBS) covering a re- 
gion of 4 GHz with a resolution of 1 . 1 MHz was used. More 
details about these observations can be found in Bujarrabal et al. 
d2010l) . The data were calibrated using the standard pipeline for 
Herschel, HIPE version 2.8. Only the H-polarization data are 
presented here because the V-polarization data are noisier espe- 
cially for the high frequency lines. We subtracted the baseline 
using a first or second order polynomial, except for the H2O 
transition at 1716 GHz, where this line is affected by standing 
waves and, consequently, the baseline was fitted using a high 
order polynomial (see Sect. 3.2). 

We detected high rotational transitions of CO as well as the 
first detection of rotational H2O in an S-star. All of the ob- 
served lines listed in Table [TJ were detected. The frequency v 
in GHz and the energy of the upper level £ up in K are given 
along with the integrated line intensity, / = jT^^dv in Kkm s -1 . 
The spectra were corrected for the main beam efficiency, i.e., 
?7mb = 0.72 exp(-(v(GHz)/6000) 2 ). The absolute calibration ac- 
curacy ranges from 10% for the lowest frequency line to 30% 
for the high frequency (> 1 THz) lines. 
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Fig. 1. The derived gas temperature structure (solid line) and the 
gas expansion velocity (dashed line) of the CSE of^ Cyg. 



3. Modelling HIFI lines 

We started the analysis by fitting the spectral energy distribu- 
tion of the CSE (assumed to be spherically symmetric) using 
Dusty (llvezic & Eli tzurlfl 9971) to derive the dust mass-loss rate. 
Based on this, we solved the gas-dust drag equation, assuming 
that both are momentum coupled, to derive the velocity struc- 
ture using the observed terminal gas velocity as a constraint, 
as shown in Fig. 13.11 This also provides a so-called dynamical 
mass-los s rate estimate, 4.9xl 0~ 7 M yr~ : in the case of x Cyg 
(see e.g., Ramstedt et al. 2008). The gas velocity law is fed into 
the CO radiative transfer model, where the line fluxes and shapes 
are computed and compared to the observations. At the same 
time, the heating by dust grains and the cooling by CO lines 
are c alculated and the resulting temperature stru cture ob tained 
(e.g., Justtanon t et al.lll994blCrosas & Ment enll997l:|Decin et al.l 
120061: Ivan der Tak et alj|2007t iRamstedt et al.ll2008l) . This pro- 
vides a circumstellar model, including a mass-loss-rate estimate 
based on the CO line modelling, which is used to model the H2O 
line emission and its contribution to the cooling. The H2O cool- 
ing is then used to recalculate the gas kinetic temperature in the 
CO model, and the process is iterated until good fits to the ob- 
served CO and H2O lines are obtained. 



3. 1 . Modelling of the CO lines 



University - Stockholm Observatory; Switzerland: ETH Zurich, 
FHNW; USA: Caltech, JPL, NHSC. 



The Monte Carlo code developed bv lSchoier & Olofssonl d2001l) 
was used to model the observed CO lines. The molecular data 
were taken from the Hitran database (Rothman et al. 2009) and 
the collisional rate coefficients from lYang et all (|2010|) for the 
41 lowest rotational levels in the v = vibrational state. We 
fitted the line shapes and fluxes for the low-/ lines obtained from 
ground-based observations, as well as the interferometric data 
for the J = 1-0 and 2-1 lines to more tightly constrain the size of 
the CO envelope (Schoier et al, in preparation). The parameters 
are listed in Table [2] The gas kinetic temperature distribution is 

shown in Fig. 13. II 

We assume a distance of 150 pc dKnapp et al.ll2003l) . and the 
resulting mass-loss rate is 7xl0 -7 M yr~' using an adopted CO 
abundance of 6x10 4 (relative to H2, see Table [2}, which pro- 
vide the best fits to both the line profile and line intensities. The 
uncertainty in the mass-loss rate is of the order 50%. This mass- 
loss rate agrees well with the dynamical mass-loss rate (certainly 
to within the errors in the input parameters). A comparison of the 
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20 

v. [km s _1 ] 




Fig. 2. The model fits (smooth lines) to the HIFI 
CO J = 6-5, 10-9 and 16-15 lines (histogram). 
The scaling factor (used to scale the model re- 
sult to fit the observed integrated intensity) is 
given for each line, showing the goodness of 
the fit. 



Table 2. Parameters used in the modelling of the line emission. 



Distance 


150 pc 


Stellar effective temperature (T c g) 


2600 K 


Stellar luminosity ( L t ) 


7.5xl0 3 L G 


Gas terminal velocity (v exp ) 


8.5 km s" 1 


Inner radius of the shell (R m ) 


2xl0 14 cm 


Gas mass-loss rate (M) 


TxlO^Meyr" 1 


CO abundance (CO/H 2 ) 


6xl0- 4 


ortho-H 2 abundance (o-H 2 0/H 2 ) 


7.5xl0- 6 


para-H 2 abundance (p-H 2 0/H 2 ) 


3.6xl0~ 6 



model fits and the HIFI observations can be seen in Fig [2] The 
models have been scaled to match the line intensities and the 
scaling factor (given in each panel of the figure) is a measure of 
the goodness of fit. The high rotational line at J — 16-15 is no- 
ticeably narrower than the lower-le vel lines observed w ith HIFI 
and ground-based instruments (e.g jKnapp et al.lll998l) . indicat- 
ing that this line originates in a region where the gas is still being 
accelerated. The estimated outer CO radius is 4xl0 16 cm. 



3.2. Modelling of the H 2 lines 

To calculate the strength and shape of the circumstellar H2O 
lines, we apply the parameters derived from our CO line mod- 
elling (mass-loss rate, gas temperature and density structure, and 
gas velocity law) to the radiative transfer mode l for H?Q based 
on an accelerated lambda iteration (ALI) code ( lJusttanont et alj 
120051: IMaercker et alj|2008l |2009|). We used the molecular data 



from JRothman et aL (|2009h and the collisional cross-sections 



from iFaure et at] d2007l) for the lowest 45 levels of ortho- and 



para-H20. The radiative excitation due to the absorption in the 
V2 bending and V3 stretching modes is included. The latter has 
been foun d to have a non -negligible effect in the low mass-loss- 
rate case ( IMaercker et al.ll2009t) . The outer ra dius of the H2O 
shell w as derived using the model results of iNetzer & Knapp 
d!987l) . i.e., 3.6xl0 15 cm. Since the ortho- and para-species are 
expected to be independent, we model the two species sepa- 
rately, using the same circumstellar input values and estimate 
the two abundances independently. Using the temperature and 
velocity structure from the CO modelling (Fig |3.1l ). we calculate 
the best-fit model to the H2O lines. All the lines are found to be 
sub-thermally excited. As mentioned above, both CO and H2O 
line cooling is taken into account. 

We present the fits to four lines of ortho-H20 and three lines 
of para-H20 observed with HIFI in Fig. [3] These lines span up- 
per energies from 60 to 300 K so the lines probe the cool part of 
the CSE as well as the warmer inner part. From Figs. l3.1l andl4l 
it can be seen that H2O lines originate well within the accelera- 
tion zone. All lines are of reasonably high signal-to-noise ratio, 
including the 3o3 -2 12 line, the highest frequency line, where the 
spectrum is affected by standing waves inside HIFI. For this line, 



we used a higher order Chebyshev polynomial for the baseline 
subtraction (bottom right panel of Fig 0). 

Our results are consistent with the velocity structure of a 
dust-driven wind as can be seen in good fits to both the CO and 
H2O line p rofiles (Figs . El and [3] respectively). Unlike the case 
for IK Tau dDecin et alj |2010a b), no modification to the dynam- 
ical calculation is required. 

4. Results and conclusions 

The observed HIFI lines are reliable probes of the inner CSE as 
the high-energy lines probe the wind in the acceleration zone. 
Both the velocity and density structures are tightly constrained 
using the HIFI lines. The abundance of CO used is 6xl0~ 4 , in- 
termediate to those usually adopted for O- and C-rich CSEs. The 
derived ortho- and para-H20 abundances are significantly lower, 
(7.5+1.4) x 10~ 6 and (3.6+0.5) x 10~ 6 , respectively (Table©. 
These values are well below the limits for O-rich AGB st ars of 
> 10~ 4 (lJusttanont et alJl2005t IMaercker et alJl2008l l2009h con- 
sistent with x Cyg being an S-star of C/O very close to unity. 
From our modelling, assuming that all carbon is locked up in 
CO (i.e., C/H = 3xl0" 4 ) and the oxygen is locked up in both CO 
and H 2 (i.e., O/H = 3x10 4 + 5.5xl0~ 6 ), our derived C/O is 
< 0.98, given that a small fractional abundance of the oxygen is 
in dust gr ains. This value is slightly higher than that of 0.95, as- 
sumed by iDuari & Hatchet! d2000l) . A non-thermal equilibrium 
chemistry model for S-stars (C/O = 0.98) predicts an H2O abun- 
dance o f 10~ 4 at the stella r photosphere, falling off to a few 10~ 6 
at 5 R* (Cherchneff 2006), an order of magnitude lower than our 
value. 

In the thermal equilibrium (TE) limit at high temperature, 
the expected ortho-to-para ratio is 3. Our derived ortho-to-para 
ratio is 2.1+0.6, close to the high-temperature TE value. The 
reported ortho-to-para ratio in CSEs of O-rich stars vary from 
1 in W H ya with a large unc ertainty dBarlow et alJI 19961) to 3 
in IK Tau dDecin et aDl2010bT) . Our result is consistent with H2O 
molecules being formed under thermal equilibrium conditions in 
the warm and dense stellar photosphere. 

Given the low total H 2 abundance of (1.1 + 0.2) x 10" 5 , 
it is clear from our analysis that the dominating cooling agent 
in the CSE of x Cyg is CO (Fig. g). Vibrationally excited H 2 
and rotationally excited H2O contribute only in the innermost 
(r < 10 15 cm) part of the CSE while CO line cooling extends fur- 
ther out until CO is photodissociated by the external UV field. 
Adiabatic cooling dominates only in the outermost part of the 
CSE. The derived H2O abundance, although much lower than in 
O-rich stars, is highe r than that obs erved in C-stars, IRC +10216 
dMelnicketalJl200ll) and V Cyg dNeufeld et alJl2010l) . consis- 
tent with AGB stars of S-type being chemically intermediate be- 
tween O-rich and C-rich AGB stars. 

Acknowledgements. HCSS / HSpot / HIPE is a joint development (are joint 
developments) by the Herschel Science Ground Segment Consortium, consist- 
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